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Abstract. Watersheds have been defined both for node and edge weighted 
graphs. We show that they are identical: for each edge (resp. node) 
weighted graph exists a node (resp. edge) weighted graph with the same 
minima and catchment basin. 



1 Introduction 

The watershed is a versatile and powerful segmentation tool. Its use for segmen- 
tation is due to Ch.Lantuejoul and S.Beucher [SJ. It may be applied on an image 
considered as a topographic surface |2], [H]> 0. An image may be considered as 
a node weighted graph ; the nodes are the pixels of the image, weighted by their 
grey tone ; the edges connect neighboring nodes and are not weighted. . Or the 
watershed may be applied on an edge weighted graph, such as a region adjacency 
graph jj, where the nodes are unweighted and represent the catchment basins, 
and the edges connecting neighboring basins are weighted by the altitude of the 
pass point separating two basins. 

In the first case, one has to find the watershed on a node weighted graphs, in 
the second on an edge weighted graph |5j . Definitions and algorithms are not the 
same in both worlds, although the have the same physical inspirations. The rain 
model, where the destiny of a drop of water falling on the surface defines the 
catchment basins ; a catchment basin is the attraction zone of a minimum, i.e. the 
set of nodes from where a drop of water may reach this minimum. Catchment 
basins generally overlap. The flooding model where the relief is flooded from 
sources placed at the regional minima and meet for forming a partition. The 
later method being often implemented as shortest distance algorithms [TO], 
There is no opposition between these models as the same trajectory may be 
followed from bottom to top, and we have a flooding model or from top to 
bottom and we have a rain model. A good review on the watershed may be 
found in [13], and in the recent book |12j . 

This paper aims at showing the equivalence between edge or node weighted 
graphs for the construction of the watershed. 



2 Graphs 



2.1 General definitions 

A non oriented graph G = [N, E] is a collection N of vertices or nodes and of 
edges E, an edge u € E being a pair of vertices (see [T],[5])- 

A chain of length n is a sequence of n edges L = {e\, &i, ■ ■ ■ , e n }, such that 
each edge e, of the sequence (2 < i < n — 1) shares one extremity with the edge 
e i-i ( e i-i 7^ 6^), and the other extremity with e^+i (ef+i ^ e$). 

A pai/i between two nodes x and ?/ is a sequence of nodes (n\ = x, n®, rik = 
?/) such that two successive nodes and n;+i are linked by an edge. 

A cycle is a chain or a path whose extremities coincide. 

A cocycle is the set of all edges with one extremity in a subset Y and the 
other in the complementary set Y. 

The subgraph spanning a set A C N is the graph Ga — [A,Ea], where E^ 
are the edges linking two nodes of A. 

The partial graph associated to the edges E' C E is G' = [N, E'] . 

A connected graph is a graph where each pair of nodes is connected by a 
path. 

2.2 Weighted graphs: regional minima and catchment basins 

In a graph G — [N, E] , edges and nodes may be weighted : ey is the weight of 
the edge (i,j) and rii the weight of the node i. The weights take their value in 
the completely ordered lattice T ■ 

Edge weighted graphs 

Regional minima A subgraph G' of an edge weighted graph G is a flat zone, if 
any two nodes of G' are connected by a chain of uniform altitude. 

A subgraph G' of a graph G is a regional minimum if G' is a flat zone and 
all edges in its cocycle have a higher altitude. 

Catchment basins A chain L = {ei, e2, ■ ■ ■ , e„} is a flooding chain, if each edge 
efc = (wfe, Ufc+i) is one of the lowest edges of its extremity riu, and if along the 
chain the weigths of the edges is never increasing. 

Definition 1. The catchment basin of a minimum m is the set of nodes linked 
by a flooding chain with a node within m. 

Node weighted graphs 

Regional minima A subgraph G' of a node weighted graph G is a flat zone, if 
any two nodes of G' are connected by a path along which all nodes have the 
same altitude. 

A subgraph G' of a graph G is a regional minimum if G' is a flat zone and 
all neighboring nodes have a higher altitude. 
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Fig. 1. On the left an edge weighted graph transformed into a flooding graph. On the 
right a node weighted graph transformed into a flooding graph. 



Catchment basins A flooding path between two nodes is a path along which the 
weigths of the nodes is never increasing. 

Definition 2. The catchment basin of a minimum m is the set of nodes linked 
by a non ascending path with a node within m. 

Remark 1. The definition of the catchment basins is the loosest possible, com- 
patible with the physical inspiration of the rain model : a drop of water falling 
on a surface cannot go upwards. With this definition, the same node may belong 
to various catchment basins. In other words there are large overlapping zones 
of the catchment basins. As most algorithm aim at producing a partition, they 
propose various methods for suppressing these overlapping zones. 

3 Outline of the method 

We want to show the equivalence of the watershed on node and edge weighted 
graphs. We first present the outline of the method on the two simple graphs in 
figjT] the left one being edge weighted and the right being node weighted. 

Consider first the edge weighted graph. It has 4 nodes A,B,C and D separated 
by weighted edges. In a flooding chain, each node is linked with the next node 
by one of its lowest edges. For this reason, if an edge is not the lowest edge of one 
of its extremities, it will never be crossed by a flooding chain. It is the case for 
the edge BC: the lowest adjacent edge of B is AB and the lowest adjacent edge 
of C is CD. For this reason, this edge BC can be suppressed from the graph 
(as presented in the second line), without modifying the flooding chains of the 
original graph. In a last step we assign weights to the nodes: each node gets the 
weight of its lowest adjacent edge, as represented in the third line .he resulting 
graph is called flooding graph. In our case, it has two regional minima, which 
are identical if one considers them from the point of view of the edge weights or 
the node weights. For this simple graph, they constitute the catchment basins. 

Consider now the node weigthed graph on the right of figJT| . It has two 
isolated regional minima. One adds a loop edge linking each isolated node with 
itself as shown in figjT] This modification does not change the flooding paths of 
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the initial graph. The last step consists in assigning, as weight to each edge, the 
maximal weight of its extremities. The added loops algo get an edge weight. As 
a result we get again a flooding graph with the following features: 

— the edges spanning the regional minima of the node weighted graph are the 
regional minima of the edge weighted graph. 

— the lowest adjacent edge of a node has the same weight as this node. For 
this reason each flooding path of the node weighted graph is simultaneously 
a flooding chain of the edge weighted graph. 

Having an identity between the minima and between flooding paths and 
chains, the catchment basins of both graphs are the same. In our case we have 
two basins with an overlapping zone containing the node F. 

4 The flooding graph 

4.1 The flooding adjunction 

We define two operators between edges and nodes : 

- an erosion [e en n] i j = tii A rij and its adjunct dilation [5 ne e] i = V e ik 

(k neighbors of i) 

- a dilation [5 en n] i j = rij V rij and its adjunct erosion [e ne e] i — A e ife 

(k neighbors of i) 

The pairs (s ne ,S en ) and (e en ,S ne ) are adjunct operators. The pairs (e ne ,S ne ) 
and (£ en ,S en ) dual operators. 

We call the first pair flooding adjunction as we may give it a physical 
explanation. Let us consider a region adjacency graph of a topographical surface, 
where rij and rij represent the flood level in the basins i and j, and e.y represents 
the altitude of the pass pont between both basins. Then: 

* the altitudes of the nodes i and j, the lowest flood covering i and j has the 
altitude [£ e n n ]jj — n i V rij 

* if i represents a catchment basin, the altitude of the pass points with the 
neighboring basin k, then the highest level of flooding without overflow through 
an adjacent edge is [e ne e]. = A e ifc 

(k neighbors of i) 

As e ne and 8 en are adjunct operators, the operator <p n = e ne 8 en is a closing 
on n and j e — S en s ne is an opening on e. 

4.2 The opening -f e 

We consider first an edge weighted graph G e and study the effect of the opening 
7 e on its edge weights. Fig{2] presents from left to right: 1) an edge weighted 
graph, , 2) the result of the erosion e ne , 3) the subsequent dilation producing 
an opening. The edges in red are those whose weight has been reduced by the 
opening. The others, invariant by the opening "f e are the edges which are, as 
we establish below, the lowest edge of one of their extremities. Two possibilities 



Fig. 2. From left to right: 1) an edge weighted graph, in the centre, 2) the result of 
the erosion e ne , 3) the subsequent dilation produces an opening. The edges in red are 
those whose weight has been reduced by the opening as they were not the lowest edge 
of one of their extremities in the initial distribution of edge weights. 

exist for an edge (i, j) with a weight A : 

* the edge has lower neighboring edges at each extremity. Hence e ne (i) < A 
and e ne (j) < A ; hence 7 e = 5 en e ne (i,j) = e en (i) V s en (j) < A : the edge is 
not invariant by the opening j e 

* the edge is the lowest edge of the extremity i. Then e ne (i) — A and 
Ene(j) < A ; hence j e = S en e ne (i,j) = e en (i) V e en (j) = A : the edge is 
invariant by the opening "f e 

Hence the edges which are invariant by 7 e are all edges which are the lowest 
edges for one of their extremities. The operator keeping for each node only 
its lowest adjacent edges is written 4- : G — > l G. As each node has at least 
one lowest neighboring edge, the resulting graph 4- G spans all the nodes. The 
resulting graph also contains all flooding chains of the initial graph (since in a 
flooding chain, each edge is the lowest edge of one of its extremities). 

The regional minima of a graph G e invariant by the opening ~f e Con- 
sider an edge weighted graph G e invariant by the opening j e . We assign to the 
nodes the weights e ne . We call G n the graph on which one only considers the 
node weights. 

Theorem 1. If an edge weighted graph G e is invariant by the opening 7 e , then 
its regional minima edges span the regional minima nodes of G n . 

Proof: A regional minimum m of the graph G e is a plateau of edges with 
altitude A, with all adjacent edges in the cocycle having a weight higher than 
A. If a node i belongs to this regional minimum, its adjacent edges have a weight 
> A but it has at least one neighboring edge with weight A : hence the weight of 
i is [e ne e) i = A. Consider now an edge (s, t), with s inside the regional minimum 
and t outside. Then e s t > X. As G is invariant by j e , the edge (s,t) is one of the 
lowest edges of the nodes t : thus the weight of t is (e ne e) t = e st > A. This shows 
that the nodes spanned by the regional minimum m form a regional minimum 
of the graph G n . 



4.3 The closing ip. 



Consider now a node weighted graph G. The closing ip n is obtained by a dilation 
S erl of the node weights followed by an erosion e ne . 

Lemma 1. The closing ip n replaces each isolated node constituting a regional 
minimum by its lowest neighboring node and leaves all other nodes unchanged. 

Proof. Consider a node with a weight A belonging to a regional minimum: 

* consider the case where the node i is an isolated regional minimum. Then S en 
assigns to all edges adjacent to i a weight bigger than A. The subsequent erosion 
e ne assigns to i the smallest of these weights, which is the weight of the smallest 
neighbor. The node i is not invariant by the closing ip n . 

* Suppose that i belongs to a regional minimum which is not isolated. The 
dilation 5 en assigns to each edge adjacent to i a weight > A. If i has a a neighbor 
j with a weight /i < A, .then d en assigns to the edge the weight A. The 
subsequent erosion £ ne assigns to i the smallest of these weights, that is A. The 
node i is invariant by the closing ip n . 

Hence if G has isolated regional minima, it is not invariant by ip n . If we add 
a loop edge linking each isolated regional minimum with itself we obtain a graph 
invariant by ip n . Indeed, if i is an isolated regional minimum with weight A, we 
add a loop edge ; the dilation S en assigns to the loop the weight A and to 
all other edges adjacent to i, a weight > A. The subsequent erosion e ne assigns 
to i the smallest of these weights, which is the weight of the loop, i.e. A. 

We write 0:G->OG the operator which adds to a node weighted graph a 
loop between each isolated regional minimum and itself. 

The regional minima of a graph G n invariant by the closing ip n Consider 
a node weighted graph G n invariant by the closing ip n . We assign to the edges the 
weights S en . We call G e the graph on which one only consider the edge weights. 

Theorem 2. If G is invariant by the closing ip n , then the edges spanning the 
regional minima nodes of G n form the regional minima edges of G e . 

Proof: A regional minimum m of G n is a plateau of pixels with altitude A, 
containing at least two nodes (there are no isolated regional minima as G n is 
invariant by ip n ). All internal edges of the plateau get the valuation A by 5 en n. 
If an edge (i, j) has the extremity i in the minimum and the extremity j outside, 
then S en n(i,j) > A. Hence, for the graph G e , the edges spanning the nodes of m 
form a regional minimum. 

4.4 The flooding graph 

We consider now a graph G on which both nodes and edges are weighted. If we 
consider only the edge weights we write G e and G n if we consider only the node 
weights. 




regional minima on the edges regional minima on the nodes 



Fig. 3. Whether one consider the egde weights or the node weights produces the same 
regional minima. 

Definition: An edge and node weighted graph G = [N, E] is a flooding 
graph iff its weight distribution (n,e) verifies both S en n = e and e ne e = n. 

In a flooding graph, the weight distribution (n, e) verifies e = 5 en n = S en £ ne e 
= 7 e e and n — £ ne e — e ne 8 en n = ip n n showing that n 6 Inv(</? n ) and e £ Inv(7 e ). 

As G is invariant by 7 e , all its edges are the lowest edge of one of their 
extremities. And as G is invariant by (p n , it has no isolated regional minimum. 

We have established earlier that: 

— if a graph G e is invariant by j e , then its regional minima edges span the 
regional minima nodes of G n . 

— if a graph G n is invariant by ip n , then the edges spanning the regional minima 
nodes of G n form the regional minima edges of G e . 

As a flooding graph is both invariant by j e and by ip n , we get the following 
theorem. 

Theorem 3. If G is a flooding graph then the node weighted graph G n and the 
edge weighted graph G e have the same regional minima subgraph. More precisely, 
the regional minima nodes of G n are spanned by the regional minima edges of 

Ge. 

Fig |3] presents the same flooding graph, on the left with its edge weights and 
on the right with its node weights : they have exactly the same regional minima. 



Flooding paths and flooding chains Consider a flooding path on G n , i.e. a 
path (rti,rt2, ■■■Wfc) of neighboring nodes with a non increasing weight, starting 
at node n\ and ending at node belonging to a regional minimum of G n . As 
Hi > Uj+x, the weight of the edge (ni,ni + i) obtained by S en is equal to the 
weight of ni and is one of teh lowest edges of r^. This shows, that the series of 
edges (ni,ni + i) form a flooding chain of G ei ending in a regional minimum of 
G e . 



Inversely, consider a flooding chain[(ni, 712), (ri2, H3), ...(rife_i, n&)] of G e . The 
weights of the edges are not increasing. Furthermore, the edge (rij, is the 

lowest edge of the node n,. As in a flooding graph the lowest edge of a node has 
the same weight than this node, the edge (n^, nj^_i) has the same weight than 
the node rij. Thus the path (ni,ri2, ■■■^fe) also is a flooding path of G„ ending a 
regional minimum of G„. 

There is a one to one correspondance between the regional minima of G e 
and G„ ; there is also a one to one correspondance between the flooding paths 
and the flooding tracks, with the same weight distribution. This shows that both 
graphs, G e on which we only consider the edge weights and G n on which we only 
consider the node weights, have the same catchment basins. 

Transforming an edge weighted graph into a flooding graph Consider 
an arbitrary edge weighted graph G e with edge weights e, and without node 
weights. The operator I G suppresses all edges which are not invariant by j e . The 
remaining edges are invariant by 7 e and verify : e = 5 en s ne e. 

If we assign to the nodes of G e weights equal to n — e ne e, then e = 5 en e ne e — 
8 en n, and the resulting graph is a flooding graph. 

Illustration Fig|4jY presents an edge weighted graph G, on which the node 
weights are produced by the erosion e ne In figj4j3 a dilation 5 en applied after the 
erosion produces an opening of the initial edge weights. The weights with a red 
color are those which are lowered by the opening. These edges are suppressed 
producing the graph J. G in fig(4p. The erosion e ne e applied on G or on I G 
produces the same node weights. The resulting graph is a flooding graph. 

Fig. [4p shows the regional minima of the complete graph if one considers 
only the node weights. In contrast, fig. [4ji! shows the regional minima of the 
flooding graph, if one only considers the node weights. They are not identical. 

If we consider only the edge weights of the flooding graph, we write G e and 
G n if we consider only the node weights. One verifies that the edges spanned by 
the regional minima nodes of G n in fig. [4^C are spanned by the regional minima 
edges of G e in figj4j^ . One also verifies that flooding paths and flooding chains 
are identical, each node being followed by an edge with the same weight. 

Transforming a node weighted graph into a flooding graph Consider an 
arbitrary node weighted graph G„with node weights n, and without edge weights. 
The operator O G adds a loop edge between each isolated regional minimum and 
itself, producing a graph invariant by (p n . The nodes verify n = e ne S en n. 

If we assign to the edges of G n weights equal to e = 8 en n, then n = e ne 6 en n = 
e ne e and the resulting graph is a flooding graph. 

Illustration Fig|5jY presents a node weighted graph. It is not invariant by tp n 
as it has isolated regional minima. One adds a loop edge linking each isolated 
regional minimum with itself, producing the graph in fig(5|3. The dilation 6 en 
produces the edge weights. The resulting edge and node weighted graph in|5p is 




Fig. 4. A: Given an edge weighted graph, the node weights are obtained by the erosioin 

B: A subsequent dilation S m produces the opening 7 e e = S en e ne e. The edges which are 
not invariant by the opening are in red. 

C: The edges which are not invariant by the opening 7 e are the edges which are not 
the lowest adjacent edge of one of their extremities. The suppression of all these edges 
produces the flooding graph. 

D: The regional minima of the complete node weighted graph do not correspond to the 
regional minima of the edge weighted graph depicted in F. 

E: The regional minima of node weighted flooding graph do correspond to the regional 

minima of the edge weighted graph depicted in F. 

F: The regional minima of the edge weighted flooding graph. 




Fig. 5. A: A node weighted graph. 

B: A loop edge is added, linking each isolated regional minimum with itself. 

C: The edge weights are produced by the dilation S en . The resulting graph is a flooding 

graph. 



a flooding graph. The regional minima have been highlighted by distinct colors. 
Again, the identity between flooding paths and flooding chains is clearly visible. 

4.5 Temporary conclusion 

We have presented how to derive from any node or edge weighted graph a flooding 
graph with the same regional minima. To each flooding path corresponds a 
flooding chain with the same weight distribution, as each node is followed by 
an edge with the same weight. The reverse also is true: to each flooding chain 
corresonds a flooding path with the same weight distribution. 

The catchment basins, as defined earlier, are exactly the same. As all non 
ascending paths and chains are accepted for defining the catchment basins, we 
also have the largest overlaping zones between them. 

Two classes of watershed algorithms have been developed, the ones for node 
weighted graphs, the others for edge weighted graphs. Since we have a one to 
one correspondance between flooding paths and flooding chains, each algorithm 
developed for a node (resp. edge) weighted graph may now also be applied for 
an edge (.resp node) weighted graph, by applying it on the associated flooding 
graph. 

5 Reducing the number of catchment basins and the size 
of their overlapping zones 

5.1 Estimating the number of catchment basins 

FigjBjAl presents a flooding graph with three regional minima. The three asso- 
ciated catchment basins overlap at a node p with weight 4. In order to break 




Fig. 6. Al : To each edge is assigned its lexicographic distance of depth 1 to the 
nearest regional minimum. A node with weight 4 has three adjacent edges with weight 
4. Keeping only one of them and suppressing the others yields 3 possible partitions 
represented in A2, B2 and C2. 

Bl: To each edge is assigned its lexicographic distance of depth 2 to the nearest regional 
minimum. The node with weight 4 has now two lowest adjacent edges with a weight 
43. Keeping only one of them and suppressing the other adjacent edges yields 2 possible 
partitions represented in B2 and C2. 

CI: To each edge is assigned its lexicographic distance of depth 3 to the nearest regional 
minimum. The node with weight 4 has only one lowest adjacent edges with a weight 
431. Keeping only this one of them and suppressing the other adjacent edges yields 1 
possible partitions represented in C2. 



the tie, only one edge adjacent to the node p should be kept, and the others 
suppressed. Three solutions are possible, illustrated by the figures [6j\2, B2 and 
C2. There is no means to decide between one or the other solution if one con- 
siders only the weight on one edge. If one considers the first two edges of each 
flooding path, we obtain a lexicographic measure, by concatenating the weight 
of the first and the second, illustrated in fig.|6j31. There remains now two choices 
between the edges with weights 43, represented in figures figures [6j32 and C2. 
Considering a third edge along the flooding paths leaves only one choice, as the 
edge with the weight 431 is the lowest. The corresponding waterhed segmenta- 
tion is presented in fig. [6p2. This example shows that we are able to reduce the 
number of partitions associated to a flooding graph, if one considers not only the 
first neighboring nodes or edges in the flooding paths or chains, but a number 
of nodes, ordered in a lexicographic order. 



By defining a lexicographic order among the flooding chains of depth 3 has 
permitted to break the ties, leaving only one solution. 



5.2 A lexicographic order relation between downwards paths 

Given a node or weighted graph we first derive its flooding graph G = [E,N]. 
We associate to G an oriented graph G by replacing the edge (p, q) by an arrow 

if n p > n q and by two arrows p% and q%> if n p = n q . The loop edge linking 
an isolated regional minimum node m with itself also is replaced by an arrow 
mrh. The graph G verifies the property (P): Vp G N, there exists at least an 
oriented path it of (with a positive or null length) linking p with a regional 
minimum. We define the catchment basin of a regional minimum as the set 
of nodes linked by an oriented path with this minimum. Obviously, each node 
belongs to at least one catchment basins. Catchment basins may overlap and 
form a watershed zone when two paths having the same node as origin reach 
two distinct regional minima. We aim at pruning the graph G, without any 
arbitrary choices, and get a partial graph G' for which the property (P) still 
holds but the watershed zones are smaller. 

As soon the path it reaches a regional minimum, it may be prolonged into a 
path of infinite length, by infinitely cycling between 2 nodes within the regional 
minimum or along the loop joining each isolated regional minimum with itself. 
All oriented paths or chains are thus of infinite length. And we may consider 
them, either in their full infinite length or consider only the first k edges. 

We now define a family of preorder relations (order relation without anti- 
symmetry) between the paths of . 

The lexicographic preorder relation of length k compares the infinite paths 
n = (PiiPi, ■■■Pki •■•) an d X = {lii ?2) ■■■Qki •••) by considering the k first nodes 
and edges: 

* 7T -< k x if n Pl < n q i or there exists t < k such that ^ < 'f ' Upi Uqi 

n vt n qt 

* 7T di k X if k ~< k X or if < k : n pi = n qi . 

This preorder relation is total, as it permits to compare all paths ; for this 
reason, among all paths linking a node p with a regional minimum, there exists 
always at least one which is the smallest for ^ fc . We say that this path is the 
steepest for the lexicographic order of depth k. 

For k = 00, we consider the infinite paths and we simply write ^. 

If 7T and% are two paths of infinite length verifying it ^ x , t nen the paths 
7T; andxi obtained by skipping the I first nodes also verify 7T; ^ xi- If 71 IS the 
smallest path linking its origin with a regional minimum, then 717 is the smallest 
path leading from pi + i to the same regional minimum. 



Nested catchment basins Consider two lexicographic order relations -< k and 
< l with I > k, then for 7Ti and 7r 2 : 7Ti ~< k 7r 2 => 7Ti -< l 7r 2 or equivalently 
h l 7T2 => 7i"2 : the steepest path for the lexicographic order I also is 



steepest for the lexicographic order k. As a consequence, a catchment basin for 
>Z l is included in the catchment basin for > k . 

For increasing values of k, the catchment basins become larger, are nested, 
and the watershed zones are reduced or vanish. For k — oo, a node is linked by 
two minimal paths with two distinct minima, only if these two paths have exactly 
the same weights, which seldom happens in natural images. In particular, if the 
regional minima have distinct weights, the catchment basins form a partition. 

Pruning the flooding graph to get steeper paths We associate to each 
order relation < k of length k a pruning operator \. k . The pruning \. k suppresses 
each edge which is not the first edge of a steepest path for < k among all paths 
with the same origin .After pruning, each node outside the regional minima is 
the origin of one or several A:— steepest flooding tracks or k— steepest flooding 
paths We say that the graph \. k G has a k-steepness or is k-steep. As for I > k, 
TTi h l 7T 2 => TTi h k 7T 2 , we have i l G C| fe G. Furthermore 4*4.' G =|'| fe G =i kyl 
G. 

Remark 2. Each pruning \. k G suppresses a number of edges still present in 
\, k ~ l G ; it suppresses them all, without doing any arbitrary choices between 
them. 

Particular k-steep graphs Applied to an arbitrary graph, the pruning | 1= | 
suppresses the edges which are not the lowest edge of one of their extremities. 
In a flooding graph, each edge it the lowest edge of one of its extremities and I 1 
is inoperant. The pruning J, 2 keeps for each node i the adjacent edges linking i 
with one of its lowest neighboring nodes. The pruning 4 00 only keeps the first 
edge of the steepest paths. 

Lemma: Any oriented path in l k G of length k is of maximal steepness for 
^ k ■ 

For this reason, a node p belongs to a k— catchment basin associated to a 
node m in a regional minimum, if there exists an oriented path in \. k G from p 
to m.For increasing values of k, the catchment basins are decreasing, and so are 
the overlapping zones between them. 

5.3 Erosions, dilations and openings on oriented graphs 

The operator | fe defined above has nice properties but is not a local operator. 
It is however possible to implement it using only local operators as we present 
now. 

Two adjunctions on oriented graphs The adjunctions (S en , e ne ) and (5 ne , e en ) 
were defined for non oriented graphs. We now define the equivalent operators 
for oriented graphs. 

The erosion from arrows to nodes assigns to each node p the minimal weight 
of all arrows having p as origin: (lt ne ) p = f\ e^. 



The dilation is obtained by adjunction. Consider a weight distribution s on 
the nodes of the oriented graph and a weight distribution e on the arrows. 

A > s p for q | p ->• q : e^ q > s p > V s p = s p = ( e n s) 

q\ P ^q q\v^q W 

The dilation (j> e nsj assigns to the arrow p% the weight of its origin p. 

We also will need the dual erosion C^enn)^ assigning to the arrow p$ the 
weight of its extremity q. 

The invariants by the opening Consider an arrow plj. The erosion s ne 
assigns to the node p the minimal weight of all arrows having p as origin. The 
subsequent dilation assigns to the arrow p% the weight of p, i.e. the the minimal 
weight of all arrows having p as origin. Thus if (7 e e )^ leaves the arrow p% 
unchanged, it means that this arrow is one of the lowest arrows having p as 
origin. 

We define a pruning operator [ ~G which cuts all arrows which are not in- 
variant by the opening ~-f e , i.e. which are not one of the lowest edges of their 
origin. 



The oriented flooding graphs We say that an edge and node weighted graph 
is an oriented flooding graph if the weights of the nodes and of the arrows verify: 

for the node weights n p = (lt ne ) p and for the edge weights: = (l) en nj _^ . 
Such a graph is invariant by the opening ~~$ en lt ne = ~-f e and by the closing 

£ ne"^ en ~^n- 



5.4 Pruning a flooding graph with local operators and without 
arbitrary choices. 

The pruning operator We start with a node and edge weighted flooding graph 
G. As explained above, we associate to G an oriented graph ~G by replacing the 
edge (p, q) by an arrow pt[ if n p > n q and by two arrows p/] and ~q% if n p — n q . 
The loop edge linking an isolated regional minimum node m with itself also is 
replaced by an arrow mrh. 

It is easy to show that if G is a flooding graph, then is an oriented flooding 
graph. 

In order to identify edges which belong to flooding chains of maximal lex- 
icographic steepness, we have to shift the weight distribution of nodes and 
weights along the oriented flooding track upwards. This is obtained thanks to 
the erosion (^ e n")^ which assigns to the arrow p% the weight of its extremity 

q- 

After this erosion, we get an edge weight distribution which is not invariant 
anymore by the opening ~rf e . Applying the pruning operator [ G leaves a graph 
which is invariant by 7 e . A final erosion (s ne ) p assigns to the nodes their new 
weights. 



Thus we have applied to ~G the following operators :~e^ en from nodes to 
arrows, the pruning operator [ on the edges followed by a last erosion e ne from 
the edges to the nodes, producing node weights n = e ne e. The resulting graph 
is an oriented flooding graph. Indeed the edges verify e = e e = ~$ en lt ne e = 
d en n. We call £ the succession of these three operators. In every day words, 
the operator £ does the following: each node is assigned the minimal weight of 
all arrows for which it is origin ; each arrow is assigned a weight equal to its 
extremity ; each arrow with a higher weight than its origin is suppressed. 

Every time that we apply £, new edges are pruned and the weight distribution 
along the flooding paths and chains moves upwards. For we obtain a graph 
which contains only flooding paths and chains of steepness > k + 1. 



Illustration 

Case of an initially edge weighted graph The fig. [7|V represents an edge weighted 
graph on which the regional minima edges are indicated in red. Fig. |7j3 presents 
the associated flooding graph G. We associate to G an oriented graph G by 
replacing the edge (p, q) by an arrow p$ if n p > n q and by two arrows p$ and q%> 
if n p = n q and get fig. [7p. The operator Q produces a new weight distribution for 
both nodes and edges in fig. [7jD. Furthermore two arrows are cut as they are not 
invariant by the opening ~^f e . As each connected component contains only one 
regional minimum, we can stop the pruning and we obtain the final watershed 
partition. This partition is indicated in false color on top of the flooding graph 
G in fig(7j5 and also in false color on top of the complete initial graph in figjTj 7 . 

Case of an initially node weighted graph Fig |8]A presents a node weighted graph. 
It contains a number of some difficulties, like the presence of 2 plateaux and of 
a buttonhole at the node with weight 3. In a first step we add loop edges linking 
each isolated regional minimum with itself as shown in fig(8j3. The next step 
produces edge weights bythe dilation S en of the node weights, yielding the initial 
flooding graph G in fLgJsJo The oriented graph G is produced in fig jsjl) . With 
each new applications of the operator £ the graph £^ is further pruned, yielding 
the graphs of fig(8]E,F and G. FigjSjT transports the arrows of the graph figjSp 
onto the graph G . The remaining flooding tracks are clearly visible, on a graph 
which has kept the initial node weights. 



6 Introducing arbitrary choices when needed 
6.1 Arbitrary choices on the flooding graph 

We now summarize the results of this study. Starting from a node or edge 
weighted graph, we want to construct a watershed partition, without arbitrary 
choices, or with a minimum number of arbitrary choices. In a first step we derive 




Fig. 7. A: An edge weighted graph. The edges in red are the regional minima edges. 
B: The associated flooding graph G 
C: The oriented flooding graph 

D: Applying the operator £ to the oriented flooding graph suppresses two edges, and 

leaves 4 connected components, containing each a regional minimum. 

E: The resulting partition superimposed on the flooding graph. 

F: The resulting partition superimposed on the initial edge weighted graph 




Fig. 8. A: A node weighted graph. 

B: A loop edge is added linking each isolated regional minimum with itself 

C: The dilation S en assigns weights to the edges and produces the flooding graph G 

C: The flooding graph is transformed into an oriented flooding graph 

E,F,G: Applying three times the operator £ to the oriented flooding graph. Each 

new application suppresses new edges, finally leaving flooding paths and chains with 

a lexicographic steepness of 4. It leaves 4 connected components, containing each a 

regional minimum. 

E,F,G: The resulting partition superimposed on the flooding graph 

F: The resulting partition superimposed with the oriented flooding graph 



the flooding graph. An edge weighted graph looses the edges which are not the 
lowest neighboring edges of one of their extremities, but keeps the edge weights. 
The node weights are obtained by the erosion e ne . To a node weighted graph 
are added some loops linking each isolated regional minimum with itself. The 
node weights keep their initial weights and the edge weights are obtained by the 
dilation S en . In both cases we obtain a graph with a perfect coupling between 
edge and node weights ; furthermore, the edge regional minima span the node 
regional minima. And the flooding chains span the flooding paths, each node 
being followed by an edge with the same weight. 

Method 1: The ties are broken by a watershed algorithm applied on the 
flooding graph. Any algorithm of the literature developed for node (resp. edge) 
weighted graphs may now be applied to the flooding graph, even if the initial 
graph is an edge (resp. node) weighted graph. The algorithm of B.Marcotegui 
et al jS] for constructing the watershed on an edge weighted graph is derived 
from Prim's algorithm for constructing a minimum spanning tree or forest. This 
algorithm is myopic and considers only flooding chains with a lexicographic 
depth equal to 1. The algorithms proposed in |5j are or the same myopic type. 
Using an algorithm [9 J based on the topographic distance [TU],[TT], and applied to 
node weighted graphs, choses flooding chains with a lexicographic depth equal 
to 2, and for this reason are more selective. Such an algorithm, designed for 
node weighted graphs may now be applied to a graph which initially was edge 
weighted. The arbitrary choices for producing a watershed partition is taken in 
charge by the watershed algorithms applied to the flooding graph. In the case 
of [5] or [5], by the scheduling of the shortest distance algorithm (ultrametric 
flooding distance for node weighted graphs or topographic distance for node 
weighted graphs). These are only a few examples of algorithms which may be 
used, among a large number of others. 

6.2 Arbitrary choices after choiceless prunings of the flooding graph 

We have seen how to reduce the number of flooding tracks and paths. The 
flooding graph G is transformed into an oriented graph ~G. The operator £ 
prunes further this graph, again without arbitrary choices. After fc iterations, 
only the flooding paths or flooding tracks of the graph G with a lexicographic 
depth fc + 1 remain. We then "transport" the arrows of the graph onto the 

flooding graph G and get a graph G^ : we keep an edge (p, q) of G if and only if 
there exists an arrow p$ or q^> in the graph £( fe ) G. Like that we obtain a flooding 
graph G in which the only remaining flooding tracks and paths have a steepness 
> k + 1. A node p belongs to two catchment basins, if there exists two flooding 
paths towards the corresponding minima, with identical fc + 1 first edges forming 
a track of fc + 1 maximal steepness. Fig(9] presents 4 pruning states of the node 
weighted graph in fig|8|\. Figj9jV presents the associated flooding graph (without 
the loops on the isolated regional minima), where the edge weights have been 
added. Two types of methods may be applied to this graph. 

Method 2: After fc steps of pruning, the graph G^ k > has only flooding paths 
and chains of steepness > fc + 1. Again any algorithm developed for node or edge 




Fig. 9. Four graphs with flooding paths of increasing steepness and decreasing over- 
lapping zones of the catchment basins. 



weighted graph may be applied to this graph. But as the remaining flooding 
paths are extremely scarce, even the loosest and most myopic algorithms will 
do a good job. In particular the algorithm by B.Marcotegui which normally 
selects flooding paths of steepness > 1 now selects paths of steepness > k + 1, 
as they are the only available. The same is true for the algorithm based on the 
topographic distance, selecting flooding paths of steepness > 2. The remaining 
arbitrary choice, if any, is taken in charge by the scheduling of the shortest 
distance algorithm (ultrametric flooding distance for node weighted graphs or 
topographic distance for node weighted graphs). 

Method 3: We continue the pruning until there remain no arrows 

which are head to tail, except in the regional minima. This is the case in figjSJ 7 . 
There exists a node with weight 4 which is the origin of two arrows with weight 
4. If we arbitrarily suppress one of them, we are also done and have a partition. 
There are two solutions possible at this stage of pruning. Thus if a node is the 
origin of several arrows, one leaves only one. (this method has often been used 
in hardware implementations, but without preliminary pruning [3], [7]). 

6.3 Maximal prunings with scarcely needed choices 

For k — oo, we obtain a maximal pruning of the flooding graph. In fact we may 
stop as soon the graph is cut into a number of components containing each one 
regional minimum. The node and edge weights remain identical as in the graph 
G, but there remains only a minimal number of flooding paths and chains. A 
node will be linked to two distinct minima by two flooding chains only if there 
exists two paths with exactly the same weight distributions towards these two 
minima ; this will rarely happen in natural images. If the minima have distinct 
weights, each node is linked with one and only one minimum by a flooding path 
of maximal steepness. 



Method 4: We continue the pruning until each connected component 

contains only one regional minimum and we are done. If this cannot be achieved, 
then we resort to method 3. After an infinite number of prunings, does 
not contain any head and tail arrows and method 3 can be applied. 

Method 5: We slightly change the weights of the minima so that they are 
all distinct. We continue the pruning until each connected component 

contains only one regional minimum and we are done. We know that this will 
happen for k < oo. 

7 Conclusion 

We have established that node and weighted graphs represent the same topogra- 
phy with the same minima and the same catchment basins. We have presented 
a method to reduce the overlapping zones of the catchment basins without arbi- 
trary choices. We finally have presented how to introduce such arbitrary choices 
if needed. 
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